Indonesian Journal of Electrical Engineering and Computer Science 
Vol. 14, No. 3, June 2019, pp. 1564~1573 
ISSN: 2502-4752, DOI: 10.1159 1/ijeecs.v14.13.pp 1564-1573 OO 1564 


Rain prediction using fuzzy rule based system 


in North-West Malaysia 


Noor Zuraidin Mohd Safar’, Azizul Azhar Ramli’*, Hirulnizam Mahdin*, David Ndzi’‘, 


Ku Muhammad Naim Ku Khalif® 


!2.3Raculty of Computer Science and Information Technology, Universiti Tun Hussein Onn Malaysia, Malaysia 
4School of Engineering and Computing, University of the West of Scotland, United Kingdom 
Faculty of Industrial Sciences and Technology, Universiti Malaysia Pahang, Malaysia 


Article Info 
Article history: 


Received Dec 15, 2018 
Revised Feb 13, 2019 
Accepted Feb 27, 2019 


Keywords: 


Expert system 

Fuzzy 

Fuzzy inference system 
Meteorology 

Rain prediction 


ABSTRACT 


The warm and humid condition is the characteristic of Malaysia tropical 
climate. Prediction of rain occurrences is important for the daily operations 
and decisions for the country that rely on agriculture needs. However 
predicting rainfall is a complex problem because it is effected by the 
dynamic nature of the tropical weather parameters of atmospheric pressure, 
temperature, humidity, dew point and wind speed. Those parameters have 
been used in this study. The rainfall prediction are compared and analyzed. 
Fuzzy Logic and Fuzzy Inference System can deal with ambiguity that often 
occurred in meteorological prediction; it can easily incorporate with expert 
knowledge and empirical study into standard mathematical. This paper will 
determine the dependability of Fuzzy Logic approach in rainfall prediction 
within the given approximation of rainfall rate, exploring the use of Fuzzy 
Logic and to develop the fuzzified model for rainfall prediction. 
The accuracy of the proposed Fuzzy Inference System model yields 72% 
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Soft computing 
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1, INTRODUCTION 

Climate change has increased incidences of unusual rainfall pattern which often cause wide spread 
flooding. Any system that can provide advanced prediction of precipitation rates could be used to predict 
flooding before the rain event. The necessity for accurate local rainfall prediction is readily ostensible when 
considering many benefits of such information would provide for river control, reservoir operations, 
forestry interests, flash flood monitoring and agriculture needs. Data fusion provides new opportunities in 
rain and flood prediction. This is vital in tropical climates where rain is localized and flash flooding is 
common. Fusion of local environmental data such as temperature, humidity, atmospheric pressure, air 
pressure, dew point, solar radiation, cloud, wind direction and wind speed from sensors could potentially be 
used to provide more accurate local predictions. 

The intention of this study is to recognize the rules and modeling rain prediction using Fuzzy Logic 
(FL) to predict rainfall. The rules will apply to forecast the rainfall in the Perlis (the northwest region of 
Malaysia). Basic environmental parameters including temperature, humidity, atmospheric pressure, 
dew point and wind speed are use in identifying the inference rules. FL is one of the core methodologies in 
soft computing which its purpose to exploit the tolerance for imprecise, uncertain, ambiguity and vagueness 
of the decision making. 
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The implementation of soft computing in engineering and science can resolve problems that 
traditionally have not been able to be solved by analytical methods. The availability of required data to make 
such predictions has been available for quite some time. The complex and changes in relationships among the 
data and its effect on the probability has often proved difficult using conventional computer analysis. The use 
of soft computing technique is to learn rather than analyses these complex relationships, has shown a great 
deal of promise in accomplishing the goal of predicting both the probability and quantity of rain in a local 
area. The methodologies in soft computing include knowledge processing and flexible knowledge acquisition 
by learning from existing data to solve non-linear problem. Study shows that FL can achieved 70% accuracy 
in data mining performance comparison [1]. The performance results are determined by assessing the 
magnitude of the error and the measurement of the correlation between the observed and forecasted values. 
Results showed that an hour ahead forecasts produce good outcomes. The outcomes indicate that Fuzzy 
Inference System models offer an ability to forecast a rainfall. This paper is structured as follows: Section 2 
describes the literature and previous study, Section 3 describes the study area and the methodology. Section 4 
presents the results and discussions and finally, Section 5 is the conclusion 


2. LITERATURE 
2.1. Related Works 

Examples on the applications of soft computing based on FL in rainfall forecasting are numerous. 
Fuzzy Rule Based System (FBS) and Fuzzy Inference System (FIS) for example have been used 
experimentally to predict weather events such as rain, fog, tornado, thunderstorm or clear condition. 
Fuzzy Time Series (FTS) presented in the time invariant and variant model for dealing with forecasting 
problems was proposed in [1]. If the annual rainfall data spanning over more than ten years, the FTS model 
can be used as an appropriate forecasting tool to predict the rainfall [2]. The design and development of the 
rules based on the three main sources; expert opinion, literature survey and automated rule generation [3]. 
The rule base is commonly obtained from an expert opinion, Handoyo and Maryji proposed the rule base 
generated based on input output pairs relationship and then each rule is optimized by ordinary least square 
(OLS) method [4]. Agboola et al. forecast the rainfall based on FIS technique. The system is entail of two 
useful modules; the knowledge base and the fuzzy reasoning unit. Two processes were executed on the FL 
model; the fuzzification and defuzzification operation [5]. Asklany et al. forecast rainfall with the help of FL 
based approximate reasoning, the IF-THEN rules. The fuzzy inference engine uses these fuzzy IF-THEN 
rules to determine a mapping from the input fuzzy sets to the output fuzzy sets based on fuzzy logic 
principles. The proposed model has the fuzzy sets contain of five environmental data: humidity, cloud, wind 
direction, temperature and air surface pressure. Those parameters are the input variables. Each fuzzy set has 
three membership functions with single output. In the fuzzy set theory, features of a fuzzy set are mapped to 
a universe of membership values using a function theoretic form belonging to the closed interval between 0 
to 1 [6]. Hasan et al. developed the rainfall prediction model using selected fuzzy variables based on the 
degree of association of different factors. The influantial in wind speed and temperature were found to have a 
positive relation with rainfall but the proposed model showed enhanced performance after adding the 
threshold values of humidity [7]. Recently, Cristani et al. use fuzzy reasoning to formulaise an architecture 
that implements the reasoning techniques in environmental study [8]. The list of rainfall prediction models 
developed using FL approached is not exhaustive. Research and achievement of the model can be found in 
[9]-[15]. The above cases are just a few examples to illustrate the potential and future direction of rainfall 
forecasting model using FL. 


2.2. Fuzzy Logic and Fuzzy Inference System 

FL is a way to definite conclusion based upon vague, ambiguous, imprecise, noise or missing input 
information [16]. FL is good at reasoning with imprecise information and explaining their behaviour based 
on fuzzy rules but they cannot automatically acquire the rules they use for making decisions. It is attempts to 
mimic human thought process. Humans analyse situations and solve problems in a relatively imprecise 
manner, do not have all the facts and the facts might be uncertain. People in a norm, generalize facts without 
having the precise data or measurements. In such situation Fuzzy approach based on FL is seems to be 
suitable. In general FL offers an inference structure that enables suitable human reasoning competencies. 
The origin of FL approach was first introduced by Lotfi Zadeh in 1965 [16]. FL is valuable in modelling 
complex and imprecise systems. Fuzzy set theory is a powerful tool and its applications have rapidly 
increased with establishing its utility in numerous areas of the scientific exploration. Fuzzy set has degree of 
membership as opposed to traditional set theory where set of membership is defined as Boolean predicate. 
Fuzzy set theory defines set membership as a possibility distribution. The general rule for this can be 
expressed as: 
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f :[0,1]' >[0,1] (1) 


Where n is some number of possibilities and n possible events and us f to generate as single possible 
outcome. The probability of FL outcome and its amount of result due to the vague input variables are being 
generated in the human mind and it is often stated as expert knowledge. This expert knowledge is the 
gathering of knowledge and ideas as a result of the experts’ experience, judgement and individual 
knowledge. In a particular system, the processes of decision making may be considered as fuzzy expressions 
seemingly adopt by the experts [7]. 

Basic structure of Fuzzy Logic model with inreference engine is presented in Figure. It is consist of 
three main steps; Fuzzification Process, Fuzzy Rules and Inference Engine Determination and 
Defuzzification Process. 
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Figure 1. Basic structure of Fuzzy Logic model 


3. STUDY AREA 

The selected area is Perlis, one of the states in Malaysia. Perlis has a 13.33 kilometre square lake 
named Timah Tasoh. There are two main rivers which are Tasoh River and Pelarit River that flows into the 
Timah Tasoh lake [17]. Agriculture area in its upstream catchments and timber plantations are surrounding 
the lake. Infrastructures development of urbanizations involved land clearing activities and amount of rain 
influenced the river water discharge of the river. Timah Tasoh Lake can be classified as medium size and 
shallow lake [18]. Existing study of environmental forecasting within stipulated area using soft computing 
approach have been conducted in [19], [20]. Three years of environmental data is recorded in hourly basis. 
The sample space for constructing the rain predictions consisted of 4624 dataset from 2012 to 2014. 80% of 
the data is used for the training and the remaining 20% for testing and validation. The training data consists 
of 153 rainfall rate reading with the minimum value of 0.02 mh"! and the maximum value of 3.10 mh”. 
Map of Malaysia and chuping as shown in Figure 2. 
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Figure 2. Map of Malaysia and chuping 
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4. APPLICATION OF RAINFALL FORECASTING 
4.1. FIS Preliminary Implementation 

FIS is designed to work with the input parameters. Once input is given by the user, appropriate rule 
for the output is given by the system. The design process for the fuzzy model will make use of the fuzzy logic 
tool box in MATLAB [21]. Figure 1s the graphical illustration of the proposed fuzzy model includes five 
input parameters (Humidity, Temperature, Pressure, Wind Speed and Dew Point) and one output parameter 
(Rainfall Rate). 
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Figure 3. Schematic illustration of the rainfall rate fuzzy model 


For each input and output variable selected, two or more membership functions (MF) is define, 
normally three but can be more qualitative category for each one of them. The proposed FIS will use four 
MFs while the output will use five MEFS. The shape of these functions can be diverse but it is usually works 
with triangles. Once the input and output variables and MEF are defined, the following step is to design the 
rule-base or decision matrix of the fuzzy knowledge-base that composed of expert IF THEN rules. 
These rules transform the input variables to an output that will tell the probability of rainfall rate. The output 
variables also have to be defined with MF, usually low, normal and high. Depending on the number of MF 
for the input and output variables the rule based will be able to define more or less potential of the existence 
of rain. The more variables existed, the more rules available to be defined, then it will make the inferences 
rules reliable. Once the realistic rules are defined according to the expert knowledge and based on collected 
weather data, these rules will become the knowledge base of each of the prediction model but it is not 
necessary to say that the whole knowledge to be translated in rules, sometimes some of the rules can be 
redundant or eliminated. This study uses the FIS Mamdani Model [22]. The Mamdani fuzzy model is similar 
to the Sugeno and Tsukomoto model where they have a set of inputs which each of the input has its own 
MFs. The Mamdani model differs from the others in the fact that the output also has a set of its own 
membership functions. In the case of the Sugeno model, the output is a function of the input MEF and in the 
case of Tsukomoto’s model; the output is a linear function of the input. 


4.2. Fuzzification 

Fuzzification is a process of transforming a real scalar value into a degree of membership (fuzzifier) 
in a particular fuzzy set value. The processes also do the mapping of the crisp input into fuzzy variable 
(linguistic label). Fuzzification of linguistic variables involves with the classification of parameters into 
linguistic labels and assigning MF for each of the variable. 


4.2.1 Determination of Membership Function 

The required data for the analysis of rainfall and construct the rules is collected from Chuping 
weather station. Three years of environmental data that consists of 4624 dataset from 2012 to 2014 have been 
used in this study. 80% of the data is used for the training and the remaining 20% for testing and validation. 
The training data consists of 153 rainfall rate reading with the minimum value of 0.02 mh"! and the maximum 
value of 3.10 mh’!. From the schematic illustration in Error! Reference source not found., five input 
parameters (Humidity, Temperature, Pressure, Wind Speed and Dew Point) and one output parameter 
(Rainfall Rate) is used to develop the fuzzy model. In order to apply Fuzzy Logic, the parameters must be 
characterized properly by flexible membership function. 


Rain prediction using fuzzy rule based system in North-West Malaysia (Noor Zuraidin Mohd Safar) 


1568 O ISSN: 2502-4752 


4.2.2 Elemntery Prediction Using IF-THEN Rules 

Elementary forms of prediction can be constructed using a policy based on simple IF THEN rules 
especially for strong inter-dependent conditions. The process of constructing the rules is to identify the 
combination of terms in the rule antecedent that best identify a class term. The combination rules can be 
implemented using the most common logical operator like AND and OR. The initial study was to identify the 
environmental parameters dependency in conjunction with rain condition. 

Cumulative Distribution Function (CDF) curves were used to determine the most significant 
parameters to predict whether it is raining or not. The CDF of each parameter is compared between rain and 
no rain events. Figureshows the curves of humidity. When it was raining, humidity varied from 68% to 
100%. Thus the following rules 1s developed, IF Humidity < 70% THEN ‘No Rain’. The rule constructed is 
IF Temperature > 33° Celsius THEN ‘No Rain’ based on temperature’s CDF. 
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Figure 4. (a) CDF for Humidity, (b) CDF for Temperature 


The following rules are constructed based on the analysis of CDF for each environmental parameter is given 
below: 


Rule A: Rule B: 
IF (Pressure >= 1014 AND IF (Pressure < 1014) OR 
Pressure < 1017) OR (Humidity > 58) OR 
(Humidity < 58) OR (Temperature < 33 ) OR 
(Temperature <33) OR (Dew Point >21 AND Dew Point <28) OR 
(Dew Point > 28) (Wind Speed > 2.5) OR 
(Wind Speed >= 4 AND THEN (Prediction event is RAIN) 


Wind Speed <=10) 
THEN (Prediction event is NO RAIN) 


The IF THEN rules prediction were constructed through analysing the training data. Rule A is used 
to predict no-rain event and Rule B to predict rain event. Rule A and Rule B are used to further develop fuzzy 
variables to feed the FIS prediction model. When rain event is detected, the ensuing question is about the rain 
intensity. In the following section the process of mapping between the fuzzy input and fuzzy output to predict 
rainfall rate is described. 


4.2.3 Input and Output Parameters 

The output parameter in Mamdani fuzzy model has a set of MF. In the proposed model a single 
output of the rainfall prediction with five MFs is expected after the appropriate fuzzification of inputs, 
inference rules and aggregation process. Details of the fuzzy model output will be discussed in the 
defuzzification process. The input parameters in this system are depicted in Table 1. The parameters is 
classified into linguistic description based on the piecewise MEFs linguistic variable representation criteria 
[23]. 
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Table 1. List of Input Parameters 


Input Parameters Unit Fuzzy Variables 
Humidity % Poor Low Medium High 
Temperature at & Poor Low Medium High 
Pressure milibar Poor Low Medium High 
Wind Speed Km/h Poor Low Medium High 
Dew Point aC Poor Low Medium High 


A fuzzy set contains elements that have varying degrees of MF in the set of ordered pairs given by: 
A=1(x,4,(x)I Vxe X)} (2) 


Where X a universal set and w(x) (usually 0 < uw,(x%) <1) is the grade of MF of the object x in A. 
A membership function 1,(«)is characterized by u,:X — (0,1) where X is the universe of discourse; x is a 
real number describing an object or its attribute and each element of X 1s mapped to a value between O and | 
[24]. MEFs can be graphically represented a fuzzy set in the form of a triangular, trapezoidal or bell shape 
(Gaussian). Every input is fuzzified by constructing the MF. In the proposed model, a triangular membership 
functions were used to describe input and output variables. Triangular membership have been used because 
of their simplicity and computational efficiency [25]. The triangular MFs are shown in Figure where a, b and 
c are the parameters of the linguistic value and the x axis 1s a value input parameter’s series. 
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Figure 5. Triangular MF 
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The MFs have been defined and built by using weather parameters and the fuzzy output is rainfall 
rate. All membership function types were selected based on statistical study of training data, synoptic and 
climitological study of the area, however most of the decision rely on knowledge and statistical analysis. 
The MFs linguistic variables for the output were adopted from the Mamdani model. In order to make the 
inference rules reliable, the rainfall rate has been classified into five linguistic labels; Poor, Low, Medium, 
High and Very High. A composition method of fuzzy relations for output has been identified using the fuzzy 
c-means (FCM) clustering method. In FCM, each data belongs to a cluster as specified by a membership 
group [26], [27]. FCM algorithm uses the reciprocal of distances to decide the cluster centres. FCM use 
weight that minimize the total weighted mean-square error. The FCM permits each feature vector to belong 
to every cluster with a fuzzy value. The algorithm assigns a feature vector to a cluster according to the 
maximum weight of the feature vector against all clusters [26]. Table and Table 2 show the MF input and 
output. Figure is a fuzzy representation based on the input output MFs for temperature. 
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Table 2. Membership Function Interval Values for Fuzzy Output Parameter 


Fuzzy Variable Rainfall Rate 
Poor 0.01 - 3.09 
Low 2.94 - 10.56 
Medium 10.05 - 2411 
High 22:93 - 46.20 
Very High 43.95 - 63.18 


Table 2. Membership Function Interval Values for Fuzzy Input Parameters 
Humidity Temperature Pressure Wind Speed Dew Point 
Poor 70.1 - 82.6 22.4 - 25.3 1001 - 1010 0.07 - 0.49 20.9 - 24.0 
Low 78.6 - 90.1 24.0 - 27.0 1005 -1012 0.47 -1.23 22.8 - 25.0 
Medium 86.3 - 96.4 25.7 - 29.4 1007 - 1014 1.17 -2.14 23.8 - 26.3 
High 91.7 - 100.6 27.9 - 30.8 1009 - 1016 2.04 - 4.41 25.0 - 26.8 
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Figure 6. Input and output variables in Fuzzy representation 


4.2.4 IF THEN Rules Development and Fuzzy Set Operations 

The most intuitive approach in developing the rules from the fuzzy set is to list all possible 
combination of the input parameter’s MF. The problems of listing all possible combination are that it makes 
the system too specific, waste system memory and increase chances of complete failure. Some rules may 
never be used. In order to reduce these problems, a proper classification technique to identify classes and 
their attribute is necessary. Decision learning tree technique is used to achieve this objective. The decision 
tree was proposed by Quinlan[28]. It is a projecting model that maps interpretations of an input to the target 
output. The induction of the rules to determine appropriate rainfall rate is an interactive process. In each 
process a new node (environmental parameter MF) is added. A node generates a number of sub-nodes until 
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all the branches of the tree belong to the same class (rainfall rate). Figure shows the illustration of the 


decision tree process beginning with the root (Humidity is High). 
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Figure 7. Fuzzy decision tree 



































4.2.5 Inferences 

According to the empirical analysis of the association between the fuzzy input and fuzzy output 
mapping and the aggregation process of the output, 37 inference rules are generated in. Some of the inference 
rules as shown in Table 4. 


Table 3. Some of the Inference Rules 


Humidity Temperature Pressure Wind Speed — Dew Point Rain Rate Prediction 

Rule 1 HIGH POOR MEDIUM POOR LOW POOR 
Rule 7 HIGH POOR MEDIUM POOR MEDIUM POOR 
Rule 11 HIGH LOW MEDIUM LOW MEDIUM POOR 
Rule 18 MEDIUM LOW MEDIUM LOW HIGH POOR 
Rule 19 HIGH POOR HIGH NONE HIGH LOW 
Rule 29 MEDIUM POOR HIGH NONE LOW LOW 
Rule 30 MEDIUM LOW HIGH NONE LOW LOW 
Rule 31 HIGH POOR HIGH LOW MEDIUM MEDIUM 
Rule 36 MEDIUM LOW HIGH MEDIUM MEDIUM HIGH 
Rule 37 MEDIUM NONE MEDIUM MEDIUM MEDIUM VERY HIGH 


4.2.6. Defuzzification 

Fuzziness benefit the system to evaluate the rules, but the final output of a fuzzy system has to be a 
crisp number [29]. It means that fuzzy numbers must be mapped to real values which can then be compared 
for magnitude. This assignment of a real value to a fuzzy number is called defuzzification. The input for the 
defuzzification process is the aggregated output from the fuzzy set meanwhile the output from the fuzzy 
system is the rainfall rate. Centroid technique is used in defuzzification process. The centre of gravity is the 
curve describing a given fuzzy quantity. Its computation requires integration of the membership functions. 
It finds the point where a vertical line would devide the aggregate set into two equal masses. Centroid 
defuzzification method finds a point representing the centre of the gravity of the fuzzy set of x, A, in the 
interval, a to b. A reasonable estimate can be obtained by calculating it over a sample of points. 
Mathematically this centre of gravity can be expressed as: 


fu, (x) xdx 
COG = ————_ (4) 


fu, (x) ax 


a 


5. PERFORMANCE INDICES 
Mean Absolute Error (MAE) and Root Mean Square Error (RMSE) were used to evaluate the 
performance of the outcome. 
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(5) 





1Y 2 
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where y, is the real value, y, is forecast value and n is the number of observations. 


The correlation coefficient between real and forecasted values of rainfall output were calculated. If 
each variable has n observations, then the Pearson correlation coefficient of R,; [30] is defined as: 


S(, ~y (>, —},) 
ee (7) 


YN n =? n = 
Ju —y,) 2 7 3° 
i=] i=l 


The R value is an indication of the relationship between the forecasted and observed rainfall. If R is 
closed to /, this indicates that there is a linear relationship and if R is close to zero, then there is non-linear 
relationship. 


6. RESULTS AND DISCUSSION 

The goal of the study is to predict rainfall rate based on the relationship between environmental 
parameters. By following the procedures described in Section 3, the procedure started with preliminary 
implementation Comparison between the actual and predicted result is important to ascertain the accuracy of 
the proposed model. The initial IF-THEN rules were used to predict rain or no-rain. The initial rules are 
important for the rainfall rate prediction. The accuracy of the prediction was evaluated using test dataset. 
The prediction for rainfall rate yielded 72% accuracy and the magnitude of error is 0.226 and 0.457 for MAE 
and RMSE respectively. 


Table shows the summarization of the results in this study. The prediction of the rainfall rate and 
actual data were carried out and compared.The results show that the FIS model is promising in predicting 
rainfall rate. It is very useful to identify the intensity of the rain using environmental parameters. 


Table 5. Results of the Performance Index 
Fuzzy Inference Model 
MAE RMSE R 
0.226 0.457 0.285 


7. CONCLUSION 

This paper introduced the intelligent model for rainfall event and its intensity prediction. 
The prediction model is rule out from the development of the FIS. Even though the finding so far not in a 
favour of reaching the appropriate accuracy for the rainfall estimation for a huge dataset, but using such 
techniques it is necessary to fuse the capabilities of expert opinion, theoretical studies and the accuracy of the 
computer systems by producing a reliable associative rules. In this model, the appropriate membership 
functions for all parameters and the process of developing efficient rules make the methodology easy to be 
driven and implemented. The models presented can be improved further by increasing the set of input 
parameters, adjusting the set of rules, hybriding with other computational methods and seasonal data. 
The results of this work suggest an open area for other environmental weather prediction such as flood, 
drought and environmental parameters. 
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